Unifying Magnons and Triplons in Stripe-Ordered Cuprate Superconductors 
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Based on a two-dimensional model of coupled two-leg spin ladders, we derive a unified picture 
of recent neutron scattering data of stripe-ordered Lai5/gBai/8Cu04, namely of the low-energy 
magnons around the superstructure satellites and of the triplon excitations at higher energies. The 
resonance peak at the antiferromagnetic wave vector Qap in the stripe-ordered phase corresponds 
to a saddle point in the dispersion of the magnetic excitations. Quantitative agreement with the 



neutron data is obtained for J = 130 — 160 meV and Jcyc/ J = 0.2 — 0.25. 
PACS numbers; 74.25.Ha, 75.40. Gb, 75.10.Jm, 75.50.Ee 
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Quantum magnetism in the cuprate superconductors 
is an intriguing issue. A detailed understanding of the 
dynamic spin susceptibility as measured by inelastic neu- 
tron scattering (INS) experiments should allow to clar- 
ify the role of magnetism in the mechanism of high- 
Tc superconductivity. In particular two features have 
been in the focus of interest: the appearance of the so- 
called resonance peak [1,2] in the superconducting phase 
at the antiferromagnetic wave vector Qaf = (1/2,1/2) 
(see Fig. lb) at finite energies (e.g. 41 meV in opti- 
mally doped YBa2Cu307_5 (YBCO)) and the existence 
of stripe order which manifests itself in superstructure 
satellites around Qaf (Fig- lb) [2-4]. In general, these 
superstructure satellites are incommensurate, but they 
may become commensurate for certain doping concen- 
trations. For many years, these two features have been 
regarded as separate issues, each of them apparent in only 
one of the two families of cuprates on which most neu- 
tron studies have focused: La2-a:Sr2,Cu04 (LSCO) and 
YBCO. But recent experimental results show that the 
resonance peak in YBCO is accompanied at lower ener- 
gies by incommensurate reflections [5,6], and that stripe 
order appears also in YBCO [7-9]. Very recently, Tran- 
quada et al. [10] observed a resonance peak at Qaf also 
in stripe-ordered Lai5/8Bai/8Cu04 (for T > Tc). 

An 5=1 collective mode (the resonance peak) in the 
superconducting phase is a prominent feature of many 
different theoretical scenarios. Its interpretation ranges 
from a particle-hole bound state (see Refs. in [1,2]) to a 
particle-particle bound state in SO (5) theory [11]. In the 
stripe-ordered phase the choice of the microscopic model 
is straightforward. Static stripe order corresponds to a 
segregation into hole-rich charge stripes and hole-poor 
spin ladders. Tranquada et al. [10] analyzed their INS 
data at high energies (including the resonance) in terms 
of the elementary triplet excitations (triplons [12]) of iso- 
lated two-leg ladders (Fig. la) which are realized in case 
of bond-centered stripes [13]. But the incommensurate 
low-energy excitations were described in a separate model 
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FIG. 1. (a) 2D square lattice split into two-leg 5=1/2 lad- 
ders (full dots) coupled by bond-centered charge stripes (open 
dots), (b) Momentum space, coordinates in reciprocal lattice 
units 27r/a (rlu). Dashed: magnetic Brillouin zone around the 
antiferromagnetic wave vector Qap = (1/2, 1/2). Open (full) 
dots: superstructure peaks from vertical (horizontal) ladders. 
(c) Momentum space rotated by 45° with h' = h + k and 
k' = k — h, thus Qap = (liO) in magnetic reciprocal lattice 
units 27r/\/2ffl (mrlu). Dotted: see caption of Fig. 4. 



as spin waves (magnons), motivated by the existence of 
weak long-range order. 

Based on a model of coupled two- leg S —1/2 ladders, 
we derive a unified description of the low-energy super- 
structure modes, of the resonance peak and of the high- 
energy excitations observed in Ref. [10]. The superstruc- 
ture modes arise from a ferromagnetic coupling between 
neighboring ladders, whereas the resonance peak corre- 
sponds to a saddle point of the triplon dispersion. The 
central result is that we arrive at a quantitative agree- 
ment with the INS data for realistic values of the ex- 
change parameters. 

The ladders shown in Fig. la are described by a stan- 
dard Heisenberg Hamiltonian (see e.g. Ref. [14]) with 
rung coupling J± > 0, leg coupling J|| > 0, and by the 
four-spin operators of a cyclic exchange 
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FIG. 2. (a), (b); dispersion ui^ and spectral weight a^(fe) 
of one trip Ion in a two- leg S—1/2 ladder with = J± = J 
and Xcyc = Jcyc/J±. (c) One-triplon dispersion LOh.k for cou- 
pled ladders at the quantum critical point J' / J — —0.072 for 

■^cvc — 0.2. 



where i counts the rungs, and R, L label the two legs. 
Inclusion of Jcyc ~ 0.2 — 0.25 is crucial in order to obtain 
quantitative agreement with experimental data, both in 
two-leg ladders [15] and in the 2D cuprates [16,17]. The 
ferromagnetic (FM) inter-ladder coupling J' < 0, see Fig. 
la, is an effective coupling which results from integrating 
out the degrees of freedom in the stripes. If two spins on 
adjacent ladders interact via a stripe rung with only one 
spin, they prefer to be antiparallel to this spin, giving rise 
to an effective FM coupling. The interaction via a half- 
filled stripe rung (two spins) implies an antiferromagnetic 
(AFM) coupling. Plausibly, the more direct FM coupling 
via one spin exceeds the AFM coupling via two spins, so 
that a weak average FM coupling results. 

We start from an effective model for isolated lad- 
ders which we have determined previously by perturba- 
tive CUTs (continuous unitary transformations) [14,18]. 
The effective model is computed perturbatively to order 
n — 10 which means that processes over a distance of up 
to 10 rungs are taken into account. The one-triplon part 
is ffl 



ladder 



n - where the operator tf (q) 
creates (annihilates) a hard-core boson - a triplon [12] - 
with spin S = 1, momentum k and one of the three triplet 
flavors a £ {x,y,z}. The dispersion w° (Fig. 2a) is ob- 
tained as a series in J\\/J± and Xcyc = Jcyc/ J i. and eval- 
uated by extrapolation [14,19]. To describe the square 
lattice we focus on Jy = where the extrapolation is 
reliable [14]. Inclusion of Jcyc reduces cj^ and in particu- 
lar the gap at A; = l/2 (by 50% for Xcyc — 0.2). 



We use the unitary transform of the observables S, 
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where the dots stand for normal-ordered quadratic and 
higher terms in the bosonic operators. On the level lin- 
ear in bosonic operators, 'S'fcg = — •S'^'S holds since rung 



triplets are odd excitations relative to rung singlets with 
respect to reflection about the centerline of the ladder. 
The Fourier transform of Eq. (2) yields 
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Below we only need the one-triplon weight a?{k) 
(Fig. 2b), which is extrapolated in the same way as uj^. 
Relative to the total weight [SfY = 1/4 of the local oper- 
ator Sf, the ^-integrated one-triplon weight a?{k) dk 
is 74% [18], 70%, 63%, and 59% for Xcyc =0, 0.1, 0.2, 
and 0.25. The missing weight belongs to multi-triplon 
channels. By inclusion of Jcyc, a^(fc) becomes more pro- 
nounced around fc=l/2. 

The key idea is to couple neighboring ladders weakly 
by J' (see Fig. la): 
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where c?/i.fe := cos{8Trh) (k) and h is the wave vector 
perpendicular to the ladders measured in rlu of the square 
lattice, cf. Fig. 1. The total Hamiltonian H is the sum of 
iJiaddcr for all ladders and H' . Since J' ^ J it is justified 
to neglect the hard-core constraint so that a Bogoliubov 



transformation yields H = J2h k-a '^h.k ^ k with 
uJh.k = \/ i^t)^ - Wdh,kOjl ■ 



(5) 



In the INS data [10] the minima are at h~\/2 ± 1/8, in 
agreement with a ferromagnetic J' < 0. 

We intend to describe the physics at the quantum crit- 
ical point where magnons just emerge [20], i.e. we choose 
J' such that ijJh,k = at the minima. The correspond- 
ing values for J' are -0.20, -0.13, -0.072, and -0.051 for 
a^cyc = 0, 0.1, 0.2, and 0.25. This is a physically plau- 
sible range [21]. The overall shape of the 2D dispersion 
LOh.k is depicted in Fig. 2c. Magnons with linear disper- 
sion emerge from (l/2±l/8, 1/2). At Qaf=(1/2, 1/2) 
a saddle point occurs. The dispersion at higher energies 
is hardly changed from its form for isolated spin ladders. 
This stems from the quadratic average (5) and from the 
suppression of the one-triplon hopping for smaller values 
of k due to the factor a^ik), cf. Fig. 2b. 

The one-triplon part of the dynamic structure factor 
Sh,k{^) reads 



^.,(.) = -^Im-^-^(-'^)«^(^)-° 



- LO 



(6) 
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where the factor a^(A:) stems from (3), the factor sin^(7r/i) 
from the interference of the two spins on a single rung, 
and oj^ from the Bogoliubov transformation. Both hor- 
izontal and vertical ladders contribute to the INS data, 
summing over different stripe domains. Thus we con- 
sider Sh^io) = (Sh^io) -|- Sk.h{^)) l"^- The constant- 
energy slices of Sh\k'{Lo) shown in Fig. 3 (in the rotated 
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FIG. 3. Constant-energy slices of Sh'^k'i^) for J = 127 
meV, Jcyc = 0.2J and J' / J = —0.072 with u ± Au given in 
each panel for direct comparison with Fig. 2 in Ref. [10]. The 
experimental Aui is accounted for hy ui ^ ui + iAui in (6). 



frame of Fig. Ic) agree very well with the experimen- 
tal results [10]. At low energies (6 meV) the cones of 
magnons emerging from the superstructure satellites de- 
fine a square around Q^f- "^he main intensity of these 
cones has shifted towards Q^p at 36 meV, culminating 
in the resonance peak at the saddle point (44 meV). At 
still higher energies, the ID character of the triplons of a 
quantum spin liquid predominates, giving rise to almost 
straight lines in the constant-energy slices. These lines 
form a diamond pattern which is rotated relative to the 
low-energy square. 

The high-energy diamond is filled in the INS data, but 
not in our single-mode calculation. This discrepancy is 
caused by the finite experimental momentum resolution, 
by the decay of a single triplon (finite life time) into two 
or three triplons due to the interladder coupling, and by 
the missing two- and three-triplon channels of the iso- 
lated ladder (see above). The two-triplon weight piles up 
around k = l/2 [18], which translates in 2D into a large 
contribution around Qaf- Energetically, the lower edge 
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FIG. 4. (a) Momentum-integrated structure factor S{u!). 
(b) S{u!) integrated only over the magnetic Brillouin zone: 
Symbols: INS data from Ref. [10]. Lines: For each value of 
a;cyc = 0; 0.1; 0.2; 0.25 with J' /J chosen to meet the quantum 
critical point (see text), we fix J such that the lower singular- 
ity peaks at 44 meV (J = 62; 85; 127; 162 meV). The intensity 
is not fitted but fixed by the value of J (using a g-factor of 
g = 2). The broadening Alj is set to 2.5meV. Thin solid lines 
(a;cyc = 0.2) mimic experimental details: (i) at high energies, 
Lorentzian broadening Au — 20meV; (ii) at low energies: in- 
tegration only in the dotted region of Fig. Ic [22]. 



of the three-triplon continuum is degenerate at Qaf with 
the one-triplon energy (e.g. for triplons with qi = Qaf, 
92 = -93 = (3/8, 1/2), and Lo{q2) =uj{q3) = 0). Thus both 
channels contribute to the filling of the diamond. 

For quantitative comparison we turn to the momen- 
tum integrated structure factor S{uj). The one-triplon 
contribution is shown in Fig. 4a. The saddle point gives 
rise to a logarithmic van Hove singularity which can be 
identified with the resonance peak observed in Ref. [10]. 
Qualitatively, pure spin wave calculations are similar in 
this respect [23]. A striking dependence of the singular- 
ity on Jcyc is found. Values of J « 140-150 meV (for 
La2Cu04) and Xcyc = 0.2-0.3 are well established [15-17]. 
Using a:cyc = 0.2 (0.25), we find J =127 (162) meV from 
the INS resonance energy. This procedure fixes also the 
absolute scale of the intensity (Fig. 4b). The quantitative 
agreement for the accepted exchange couplings strongly 
corroborates our approach. The remaining differences 
can be explained easily, (i) Our theory predicts a second 
sharp peak above 200 meV. There, the experimental res- 
olution is only ±20 — 25 meV, so that the high-energy 
peak is washed out entirely (thin solid line in Fig. 4b). 
(ii) We find too much weight below « 40 meV because 
we do not consider a finite staggered magnetization. 

The superstructure periodicity of 4a with antiphase 
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domain structure manifests itself in the reciprocal space 
coordinates of the satellites (n/8, 1/2) and (1/2, n/8) 
with n odd. This order has been depicted either as three- 
leg ladders separated by a sitc-ccntcrcd charge stripe or 
as the bond-centered charge stripes discussed here and in 
Ref. [10]. Recent ab initio results [13] support the bond- 
centered scenario. The qiiantitativc agreement in Figs. 3 
and 4 corroborates the two-leg ladder scenario, but cou- 
pled three-leg ladders cannot be ruled out at present. 

Earlier, a unified description of the resonance peak and 
of the low-energy satellites was proposed based on strictly 
ID models [24]. But in ID the satellites carry much more 
spectral weight than the resonance. A true resonant be- 
havior emerges only from a 2D saddle point. 

In stripe-ordered 5 = 1 La2-xSr^Ni04 (LNO) [25,26], 
isotropic magnons emerging from the superstructure 
satellites describe the INS data up to the highest mea- 
sured energy (w 100 meV) [23]. In the nickelates, a spin 
wave picture is much more justified due to the larger spin 
S=l and due to the large inter ladder coupling J', which 
was found to be J' /J « 0.5 - 1 [25,26,23]. This shifts the 
saddle point, where the magnon cones merge, to high en- 
ergies. Indeed, enhanced intensity was observed at ~ 80 
meV at Qaf [25]. In the cuprates, the saddle point is 
rather low due to the small value of J' and due to the 
sizeable value of Jcyc (Fig. 4a). 

Our spin-only model aims at the magnetic response of 
the stripe-ordered phase. Let us nevertheless speculate 
about the possible implications for the magnetic proper- 
ties of the superconducting phase. Since our calculation 
does not rely on long-range magnetic order, it is plausible 
that its qualitative features survive in the superconduct- 
ing phase. The phase on the disordered side of the quan- 
tum critical point shows a finite spin gap [20] . as observed 
experimentally [5,6,9]. A saddle point exists if the dis- 
persion is large in one direction and small in the perpen- 
dicular direction. But in a metallic phase the damping 
due to charge degrees of freedom has to be considered. In 
stripe-ordered Lai5/8Ba]^/8Cu04 the spin-only model is 
applicable because the damping is reduced by the charge 
order. Similarly, the appearance of the peak and of the 
incommensurate low-energy modes below in nearly op- 
timally doped YBCO [5] can be attributed to the van- 
ishing damping if the resonance mode is lying below the 
particle-hole continuum. Note the stunning similarity of 
Fig. 3 with recent INS data of YBCO [27] showing de- 
pleted diamonds. The absence of the resonance peak in 
LSCO, with similar J but a smaller superconducting gap, 
suggests that the mode is located within the continuum, 
where it is overdamped [28,29]. The downward disper- 
sion below the resonance energy iVr [5-7,9] as well as the 
upward dispersion above iVr [6,7,9] are generic features of 
our stripe model. Close to ojr the saddle point implies an 
almost circular intensity distribution in Sh,k which may 
reconcile recent INS data of optimally doped YBCO [6] 
with a stripe scenario. 



In conclusion, we have given a unified description of 
low-energy 2D magnons, high-energy ID triplons and 
of the resonance mode at Qaf in stripe phases. Both 
the resonance energy and the spectral intensity agree 
quantitatively with the INS data of Lai5/8Bai/8Cu04 
[10] for realistic values of J w 130 — 160 meV and 
Jcyc/ J ~ 0.2 — 0.25. This imderlines the significance 
of quasi-lD quantum magnetism in the 2D cuprates. 

We thank J.M. Tranquada for provision of the INS data 
and acknowledge fruitful discussions with M. Braden and 
the support by the DFG via SFB 608 and SP 1073. 

While we wrote this manuscript we became aware of 
a mean-field calculation of cu^.k based on coupled dimers 
[30] leading to qualitatively similar results as in Fig. 3. 
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